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Abstract 

Context. The height structure of 3-min oscillations over sunspots is 
studied in the context of the recently discovered effect of height inversion: 
over the umbra, the spatial location of the maximum of chromospheric 3- 
min oscillation power corresponds to the relative decrease in the power of 
photospheric oscillations. 

Aims. We investigate whether the height inversion of the power of 3-min 
oscillations is a common feature of the spatial structure of the oscillations 
for the majority of sunspots. 

Methods. Spectrogram sequences of Ha 6563A and Fe i 6569A over 
sunspots, acquired with very high cadency (about 2 s or better) are obtained. 
The distribution of the oscillation power of the line-of-sight velocity sig- 
nal is studied by using methods of wavelet frequency filtration and Fourier 
analysis. 

Results. The efi'ect of the height inversion is found in 9 of 1 1 analyzed 

active regions. The interpretation of this effect is possibly connected to both 
the decrease in the level of photosphere in sunspot umbrae and the magnetic 
field topology. 
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1 Introduction 



Powerful 3-min oscillations of the line-of-sight (LOS) velocity, detected in the 
Doppler shifts of the chromospheric spectral lines above sunspots have been sub- 
ject to intensive studies for several decades (Beckers and Tallant, 1969; Giovanelli, 
1972; Zirin and Stein, 1972; Lites, 1992). The interest in this phenomenon is jus- 
tified for several important reasons. First of all, it is connected to the transfer of 
the energy of photospheric motions by these waves to the corona. There is obser- 
vational evidence of 3-min field-aligned compressible waves in the corona over 
sunspots, usually seen in 171 A and 195A passbands of TRACE (see, e.g. (De 
Moortel, 2006)), which may be associated with the leakage of the 3-min sunspot 
oscillations. As the waves are interpreted as slow magnetoacoustic waves that 
propagate almost parallel to the magnetic field, their path highlights the mag- 
netic connectivity of the photospheric and coronal plasmas. In addition, these 
waves are a useful tool for probing solar coronal plasmas in the quickly devel- 
oping research field of coronal seismology (King et ah, 2003; Wang, Innes, and 
Qiu, 2007). Moreover, the relationship between 3-min oscillations in sunspots 
and the appearance of a similar periodicity in flaring energy releases has been es- 
tablished (Sych et ah, 2009). Another useful application of 3-min oscillations is 
in attempting to reveal the internal, sub-photospheric structuring of sunspots (e.g. 
(Zhugzhda, 2008)). 

3-min sunspot oscillations are presently studied in the optical band by both 
filtering and spectral methods, as well as in radio and UV bands (Gelfreikh et al, 
1999; O'Shea, Muglach, and Fleck, 2002; Rouppe van der Voort et al, 2003; 
Bloomfield, Lagg, and Solanki, 2007; Nagashima et al, 2007; Jess et al, 2007; 
Tziotziou etal, 2007; Balthasar and Schleicher, 2008). (Kobanov and Makarchik, 
2004; Kobanov, Kolobov, and Makarchik, 2006; Bloomfield, Lagg, and Solanki, 
2007) demonstrated that 3-min oscillations in the umbral chromosphere are not 
the source of the running penumbral waves and concluded that the apparently 
horizontally propagating waves correspond to the "visual pattern" scenario. 

Another identified peculiarity in the behaviour of 3-min oscillations has been 
the effect of height inversion (Kobanov, Kolobov, and Chupin, 2008): at the chro- 
mospheric level, the spatial localisation of the maximum of 3-min oscillation 
power, usually over the umbra, often corresponds to the minimum of these os- 
cillations at the photospheric level. One would expect that the spatial distribution 
of 3-min oscillations at both heights is similar, although at the photospheric level 
the oscillation power in the umbra is lower than in the penumbra and adjacent re- 
gions. This finding is consistent with the results of previous studies. In particular. 
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it was pointed out that in the umbra, at the photospheric height, 3-min oscillations 
are either not detected at all (Balthasar, Kueveler, and Wiehr, 1987) or of very low 
amplitude relative to the level of noise (Lites and Thomas, 1985). We note that 
the suppression of the photospheric oscillation power in active regions has also 
been detected by the helioseismological methods (Braun and Duvall, 1990; Braun, 
1995; Nicholas, Thompson, and Rajaguru, 2004). The analysis of the brightness 
fluctuations observed in the G-band with Hinode/SOT detected the decrease in the 
broadband oscillation power in central parts of sunspots (Nagashima et al, 2007). 

The aim of this paper is to perform a comparative analysis of 3-min oscilla- 
tions for a number of sunspots, demonstrate that the height inversion is a statisti- 
cally significant and reproducible effect, and contribute to its interpretation. Our 
approach is based on the use of high cadence optical data, with the time resolu- 
tion of about 2 s, and upon the consideration of 3-min oscillations in a narrow 
band, 4.7-6.7 mHz. The paper is organised as follows. In the next section, we 
describe the instrument used in the observations and data reduction. In Section 
3, we present the results obtained. The results obtained are summarised in the 
conclusions. 

2 The instrument and data reduction 

We use data obtained with the Horizontal Solar Telescope at the Sayan Solar Ob- 
servatory, Russia (Kobanov et al, 2009). The diameter of the mirror of the instru- 
ment is 80 cm, hence the theoretically possible spatial resolution is 0.2". However, 
because of the Earth's atmospheric conditions the resolution is usually about 1". 
The guiding system carries out targeting and object capturing with a precision no 
worse than 1". In observations, an astronomical CCD of the resolution 256x1024 
pixels was used. The spectral snapshot contains information about a spatial re- 
gion of the size 63x0.5", defined by the slit of the spectrograph. For each spatial 
element, we obtain the spectrum of the width of 8. 5 A (for the spectrograph dis- 
persion of the fifth order in the vicinity of the Ha 6563 A line). The image was 
rotated in a way to position the spectrograph slit at the sunspot centre in either 
the north-south or east- west directions (see Figure 1), by using a Dove prism in 
front of the spectrograph. The duration of observational sequences was about one 
hour with a cadence time of about two seconds. The majority of the observational 
sequences were recorded with the use of a deflector (Kobanov, 2001), which pro- 
vided us with information about the LOS velocity, spectral line intensity, and the 
LOS component of the photospheric magnetic field (by the Fe i 6569A line) si- 
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multaneously. 

The information was recorded in a digital data cube of the following structure. 
Two dimensions corresponded to the pixel number along the spectrograph slit and 
dispersion (hence the spatial and spectral information), while the third dimension 
contained the intensity of the emission at the certain instant of time. A dedicated 
software package was then applied to reduce the data cube to a two-dimensional 
array. The array contained the distribution of the intensity and the LOS compo- 
nents of the velocity and magnetic field along the slit at each instant of time. In 
the analysis of oscillatory processes, the time sequences of each spatial pixel were 
studied by the Fourier and wavelet techniques. The wavelet analysis was per- 
formed according to the methodology described in Torrence and Compo (1998) 
with the 6th-order Morlet basis function. 

The LOS component of the velocity was obtained by the lambdameter tech- 
nique (Rayrole, 1967). In this method, the spectral position of a chosen line is 
determined by two virtual slits. The distance between the slits is known and re- 
mains constant during the measurements. Initially, the slits are situated at equal 
distances from the line centre, hence have equal intensities. If in the next spectro- 
gram the line is displaced, the intensities measured by the slits differ. Displacing 
the slits to make their intensities equal allows one to determine the new location 
of the line. Since this method is based upon the determination of the relative posi- 
tion of the line, it is not affected by the variation in the intensity in time, from one 
spectrogram to another. 

In the processing of the observational sequences obtained with the deflector, 
the LOS velocity is determined by the Doppler formula directly from the displace- 
ment of He line, as the splitting of the line by the deflector is insignificant. In the 
case of another line used in this study, Fe i, it is completely split by the deflector 
into the cr_ and (t+ components. The LOS velocity is proportional to the sum of 
the displacements of these components, while the magnetic field is proportional 
to the difference in the displacements. The magnetic field is then calculated ac- 
cording to the Zeeman effect for a simple triplet. Telluric line near Ha was used 
to eliminate the spectrograph noise. 

A possible error in the determination of the line position may be related to the 
asymmetry of the line shape. The asymmetry can be caused by blending of the 
He line, or by the variation in the LOS speeds between the heights covered by the 
optically thick Ha line. However, as in this study we consider not the absolute 
values of the LOS speeds, but the frequencies of their time variations, the latter 
effect not being likely to cause significant error. In this paper, we did not discuss 
the results of the magnetic field measurements. 
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Figure 1 : An example of the spectrograph slit orientation in the observations of 
the sunspot NOAA 791 (the scale of the slit is not kept). 




Figure 2: The distribution of the wavelet power of the LOS velocity in the range 
4.7-6.7 mHz in the chromosphere of NOAA 810 observed in Ha. Left panel: 
time-spatial dynamics of the 3-min oscillations. Central panel: Profile of the time- 
averaged spatial distribution of the wavelet power. Right panel: the profile of the 
intensity observed in the continuum. The vertical dashed hnes indicate the umbra- 
penumbra boundary. 

3 Results 

In our study, we analysed 1 1 active regions (AO) observed during 2002-2007. A 
list of the AO with their locations on the solar disk, the time of the observation, 
and results is presented in Table 1 . The duration of each sequence is about one 
hour, where Pp and Pc are the ratios of the mean power of the 3 min oscillations 
in the spot umbra to the mean power in the penumbra for the photospheric (Pp) 
and chromospheric (Pc) level correspondingly. 

The spatial localisation of the 3-min oscillations of the LOS velocity has been 
studied using both Fourier and wavelet analyses. In both cases, the analysed spec- 
tral interval was taken to be from 4.7 to 6.7 mHz. Figure 2 shows a typical wavelet 
spectrum that gives the time-spatial dynamics of the 3-min oscillation, as well as 
the time-averaged spatial distribution of the oscillation power. 
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Table 1 : Observational sequences studied in this paper. 



NOAA 


Solar coordinates 


Date 


Start time, UT 


Pp/Pc 


051 


S17E01 


02.08.2002 


00:45 


0.5/12 


051 


S16E10 


01.08.2002 


06:52 


0.8/12.9 


105 L 


S08E39 


11.09.2002 


01:21 


0.6/9.2 


613 


S09W03 


20.05.2004 


03:25 


0.6/12.2 


621 


S14W06 


04.06.2004 


05:43 


— 


657 


N10W03 


13.08.2004 


01:44 


0.6/7.3 


661 L 


N06E42 


15.08.2004 


23:57 


0.7/3.7 


661 L 


N07E20 


18.08.2004 


00:16 


0.7/11.4 


661 F 


N07E19 


18.08.2004 


01:01 


1.2/38.2 


791 L 


N13E21 


26.07.2005 


01:37 


0.5/5.1 


791 L 


N13E20 


26.07.2005 


03:27 


2.1/15.6 


791 L 


N13E08 


27.07.2005 


00:53 


0.3/26.5 


791 L 


N13E06 


27.07.2005 


01:26 


0.5/16.3 


791 L 


N13E06 


27.07.2005 


03:06 


0.6/25.3 


791 L 


N13E05 


27.07.2005 


04:46 


0.4/33.3 


794 


S11E32 


04.08.2005 


02:47 


0.7/11 


794 


S11E32 


04.08.2005 


03:38 


1.1/18.6 


794 


S11E05 


05.08.2005 


01:22 


0.8/93 


810 


N10E23 


21.09.2005 


02:45 


0.9/9.2 


810 


N10E23 


21.09.2005 


06:48 


0.8/6.3 


810 


N09W07 


23.09.2005 


07:15 


1.5/9.2 


886 L 


N07E24 


25.05.2006 


00:47 


0.8/3.5 


886 L 


N08E02 


26.05.2006 


04:25 


0.7/4.7 


963 L 


S07E62 


08.07.2007 


23:37 


0.5/4 
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Figure 3: The effect of the height inversion for the power of 3-min oscillations in 
sunspot umbrae. The panels show time- averaged spatial profiles of 3-min oscilla- 
tions of the LOS velocity measured in different active regions. Thin lines corre- 
spond to the photospheric signal (Fe i), and the thick curves to the chromospheric 
signal (H a). The vertical dashed lines indicate the umbra-penumbra boundary. 
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Figure 4: The absence of the height inversion effect for the leading sunspot of 
NOAA 791 in the observational sequences on 2005 July 26, 03:27-5:40 UT. In 
the same sunspot, the effect of the height inversion is presented in other five time 
series (one of them is shown in Figure 3). The figure shows the spatial distri- 
bution of the wavelet power of the LOS oscillations in the band 4.7-6.7 mHz. 
Top panels: Time-spatial dynamics of the wavelet spectra. Bottom panels: Time- 
averaged profiles of the wavelet spectra. The vertical dashed lines indicate the 
umbra-penumbra boundary. 
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Our analysis demonstrated that at the chromo spheric level in all considered 
sunspots except the small sunspot NOAA621, the maxima of 3-min oscillation 
power are situated near the centres of the umbrae. At the photospheric level, in 
the most considered time series, these regions were found to correspond to the 
suppression of the 3-min oscillation power (see Figure 3). The criterion used in 
the detection of the effect of the height inversion is the occurrence of a local min- 
imum in the spatial distribution of the photospheric narrowband power and its 
absence at the chromospheric level. The spatial coordinates in all plots of Fig- 
ure 3 were restricted by the external boundaries of the penumbra. The effect can 
be quantified by the ratio of the mean power in the umbra (averaged over the whole 
umbra) to the mean power in the penumbra (averaged over the whole penumbra). 
In the chromosphere, this ratio was found to be greater than one, while in the pho- 
tosphere it is smaller than one (see Pp/Pc, Table 1). We analysed 24 time series 
from 1 1 active regions. In 18 series, we found the effect of height inversion, in 3 
series the situation was the opposite, and for 2 series we derived uncertain results 
(the photospheric mean power ratio was 0.9-1.1). One time series of NOAA 621 
shows the absence of the 3-min LOS velocity oscillations in the chromosphere. 
We should also point out that in certain time intervals in some sunspots the maxi- 
mum of 3-min oscillations at the photospheric level can be found in the umbra. For 
example, for NOAA 791, in six analysed time sequences, the effect of the height 
inversion was found in five cases, while in one case the photospheric maximum of 
the 3-min oscillations was localised in the umbra, hence exactly coincident with 
the chromospheric maximum (see Figure 4). Wavelet spectra of the photospheric 
and chromospheric signals are shown in the upper panels. At first glance, the 
increases and decreases in the oscillation power at those heights seem to occur 
simultaneously. However, the detailed phase analysis discussed below does not 
support this impression. 

We are thus inclined to conclude that in the majority of the observed cases the 
spatial location of the maxima of 3-min chromospheric oscillations correspond to 
the suppression of these oscillations at the photospheric level. One possible inter- 
pretation of this effect may be the decrease in the photospheric level observed in 
Fe I 6569A (the Wilson depression), as it is known that the amplitude of acoustic 
oscillations decreases in deeper regions (Lites, 1992; Sigwarth and Mattig, 1997). 
The question is whether the difference in heights, associated with the Wilson de- 
pression, is responsible for the observed decrease in the oscillation amplitude by 
a factor of two. 

According to Balthasar and Woehl (1983) and Watson et al. (2009), the depres- 
sion (t=1) of the umbral photosphere with respect to the penumbra (the Wilson 
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effect) is about 700 km. This difference in height can lead to the decrease in the 
3-min oscillation amplitude by a factor of 2 (Lites et ah, 1998). This decrease 
is consistent with the 3-min LOS velocity power distribution at the photospheric 
level (Figure 3), at least within the order of the magnitude. 

Pamell and Beckers (1969) used the formation depth of the line Fe I 6569 
to study the granulation. They determined the upper limit to the line formation 
height to be 250 km. According to the broadly used atlas of the solar spectrum 
(Moore, Minnaert, and Houtgast, 1966), this line is suppressed in sunspots, hence 
one can assume that its optical thickness in the umbra is not higher. In the Fe I 
6569 line, the Wilson effect is as pronounced as in the continuum. 

The observed effect may also be caused by the topology of the sunspot mag- 
netic field. If we assume that a horizontally extended source of 3-min oscillations 
is situated below the photosphere, in the umbra, the oscillations in the form of slow 
magnetoacoustic waves would be guided by the vertical magnetic field lines up- 
wards and reach the chromosphere. At the edge of the umbra and in the penumbra, 
the magnetic field lines are inclined from the vertical direction, preventing the free 
propagation of the short period waves (Bogdan and Judge, 2006). Thus, under the 
canopy 3-min waves are able to propagate in the vertical direction up to the mag- 
netic dome only. Finsterle et al. (2004) pointed out that the 7-mHz oscillations are 
reflected by the active region canopy. The partial reflection of the waves from the 
magnetic dome can lead to the observed accumulation of the wave energy at the 
photospheric level under the regions with a horizontal magnetic field. This may 
explain the relative increase in the 3-min photospheric oscillation power in the 
regions outside the umbra, and the corresponding decrease in the chromospheric 
oscillation power in these regions. 

Moreover, it is necessary to consider a possible contribution of the "power 
halo"(Lindsey and Braun, 1999; Braun and Lindsey, 1999). Our observations 
were focused on the sunspots, and, unfortunately, studying the sunspots' vicinities 
is beyond the scope of the data for most time series. 

In the majority of the observational sequences, the power of photospheric 3- 
min oscillations in the umbra was low, thus it was impossible to estimate the 
phase delay between the photospheric and chromospheric oscillatory motions in 
the umbra. However, the above-mentioned unusual observational sequences of 
NO A A 791, containing the umbral maximum of the photospheric 3-min oscilla- 
tions, allows us to carry out the analysis of the phase shifts between the signals 
at different heights. Unfortunately, the direct comparison of the photospheric and 
chromospheric signals does not provide us with any conclusive results. Figure 5 
shows that there is no stable phase difference between narrowband photospheric 
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time, min 

Figure 5: Left panel. Comparison of phases of chromo spheric (thick line) and 
photospheric (thin line) narrowband (4.7-6.7 mHz) LOS velocity signals in the 
centre of the umbra of NOAA 791 (the leading spot) in different time intervals. 
Right panel. Correlation between the chromospheric and photospheric signals of 
the LOS velocity for the whole time series. 

and chromospheric signals. In some time intervals (e.g. 45-50 min), the chromo- 
spheric signal even precedes the photospheric signal. There are also time intervals 
(e.g. 1 12-130 min) when the time lag gradually increases from 20 s to 140 s. The 
correlation between the signals for the whole time series is presented in the right 
panel of Figure 5. 

A similar negative result comes from the comparison of the power maxima 
of the photospheric and chromospheric signals (Figure 6). The detuning of the 
signals is so pronounced that one has the impression that the photospheric and 
chromospheric signals are disconnected. Thus, the data obtained does not allow 
us to measure the phase lag between the photospheric and chromospheric 3-min 
oscillations directly, to derive the phase speed. However, the average time lag be- 
tween the appearance of the maxima of the oscillation power at those two heights 
could be measured and is about 150 s. Taking the distance between the heights 
of observations to be 2000 km (Vernazza, Avrett, and Loeser, 1981; White and 
Wilson, 1966), we obtain the vertical group speed of about 13kms"^ 
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Figure 6: The time variation in the power of the chromo spheric (thick line) and 
photospheric (thin line) narrowband (4.7-6.7 mHz) LOS velocity signals in the 
umbra of the leading spot of NOAA 791. The upper panel shows the signals of 
the spatial element of 1" at the centre of the umbra. The bottom panel shows 
the signals averaged over the 10" region in the umbra. Right panel shows the 
correlation between the signals presented in the upper panel. 

On the other hand, taking into account the correlogram shown in the right 
panel of Figure 6, the group speed can be estimated as 6km/s. This discrepancy 
is due to the ambiguity of the phase connection between the photospheric and 
chromospheric signals. Moreover, the observed correlation coefficient is rather 
low, about 0.35, and is observed only in one case. Thus, the present study does 
not allow us to make confident and statistically- significant conclusions about the 
wave connectivity at the photospheric and chromospheric levels, based upon the 
correlation analysis. One possible reason for the apparent out-of-phase behaviour 
of the photospheric and chromospheric oscillations may be the partial reflection 
of the slow magnetoacoustic waves from the density gradient over the umbra, 
which causes a complicated time-dependent superposition of the propagating and 
standing waves in the photosphere. 
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4 Conclusions 



We have analysed the spatial distribution of narrowband 3-min oscillations of the 
LOS velocity at the photospheric and chromo spheric levels in eleven sunspots. 
The observations dedicated to this study were carried out with the 80 cm Hor- 
izontal Solar Telescope at the Sayan Solar Observatory in the Fei 6569A and 
Ha 6563A lines, respectively. The spectral resolution of the instrument is 8.5 mA. 
The cadence time of the measurements was about 2 s. At the chromospheric level 
in ten of the eleven considered sunspots, the maxima of 3-min oscillation power 
were found to be situated near the centres of the sunspot umbrae. In nine of the 
eleven sunspots at the photospheric level, the 3-min oscillation power was found 
to be suppressed in the umbra relative to the peripheral regions. 

Thus, in the umbra, the chromospheric oscillations have higher power than 
outside, while, in contrast, the photospheric oscillations in the umbra have sys- 
tematically lower amplitude than in the surrounding regions. This effect may be 
termed a height inversion of 3-min oscillations in sunspots. Phase shifts between 
narrowband 3-min amplitude photospheric and chromospheric signals were found 
to vary in time, as well as between the maxima of their power. The effect of height 
inversion may be attributed to the wave reflection caused by the inclined mag- 
netic fields in the peripheral regions of sunspots or to the Wilson depression of 
the photospheric level in the umbra. However, detailed theoretical study of these 
phenomena is needed. 
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